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Abstract 

Background: Blood-brain barrier (BBB) disruption is an integral feature of numerous neurological disorders. However, 
there is a relative lack of knowledge regarding the underlying molecular mechanisms of immune-mediated BBB 
disruption. We have previously shown that CDS T cells and perforin play critical roles in initiating altered permeability of 
the BBB in the peptide-induced fatal syndrome (PIFS) model developed by our laboratory. Additionally, despite having 
indistinguishable CDS T cell responses, C57BL/6J (B6) mice are highly susceptible to PIFS, exhibiting functional motor 
deficits, increased astrocyte activation, and severe CNS vascular permeability, while 129Sl/SvlmJ (129S1) mice remain 
resistant. Therefore, to investigate the potential role of genetic factors, we performed a comprehensive genetic analysis 
of (B6 X 1 29S1) F2 progeny to define quantitative trait loci (QTL) linked to the phenotypic characteristics stated above 
that mediate CDS T cell-initiated BBB disruption. 

Results: Using single nucleotide polymorphism (SNP) markers and a 95% confidence interval, we identified one QTL 
[PlFSl] on chromosome 12 linked to deficits in motor function (SNP markers rs6292954, rsl34S1303, rs3655057, and 
rsl34S1324, LOD score = 3.3). In addition we identified a second QTL (P/F52) on chromosome 17 linked to changes in 
CNS vascular permeability (SNP markers rs6196216 and rs3672065, LOD score = 3.7). 

Conclusions: The QTL critical intervals discovered have allowed for compilation of a list of candidate genes implicated 
in regulating functional deficit and CNS vascular permeability. These genes encode for factors that may be potential 
targets for therapeutic approaches to treat disorders characterized by CDS T cell-mediated BBB disruption. 

Keywords: Quantitative trait loci (QTL), Single nucleotide polymorphism (SNP), Blood-brain barrier (BBB), Theiler's 
murine encephalomyelitis virus (TMEV), Peptide-induced fatal syndrome (PIFS), CDS T cell, CNS vascular permeability 



Background 

Disruption of the blood brain barrier (BBB) is a common 
feature of several devastating immune-mediated neuro- 
logical disorders, including multiple sclerosis, dengue 
hemorrhagic fever, cerebral malaria, stroke, and acute 
hemorrhagic leukoencephalitis (AHLE) [1-8]. The under- 
lying molecular mechanisms of BBB disruption observed in 
these conditions remain poorly understood, undermining 
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the development of potential therapeutic approaches for 
human diseases characterized by BBB disruption. Our lab 
developed an inducible murine model of CNS vascular 
permeability using a variation of the Theiler s murine en- 
cephalomyelitis virus (TMEV) model commonly used to 
study multiple sclerosis. In this model, peptide-specific 
stimulation of CNS -infiltrating CDS T cells in vivo initi- 
ates disruption of the BBB, resulting in astrocyte activa- 
tion, severe CNS vascular permeability, and functional 
motor deficits, followed by death within 48 hours. This 
peptide induced fatal syndrome (PIFS) is dependent on 
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perforin expression, but does not rely on neutrophils, CD4 
T cells, TNF-a, IFN-y, LT|3R or IL-1 [9,10]. 

In addition to CNS- infiltrating CDS T cells and their 
effector functions, there is potential for other genetic 
factors to contribute to BBB disruption. C57BL/6J (B6) 
and 129Sl/SvImJ (129S1) mice both express the H2-D^ 
class I molecule and thus are resistant to chronic TMEV 
infection and demyelinating syndrome [11-13]. However, 
despite having identical MHC class I molecules, epitope- 
specific CNS -infiltrating CDS T cells, and equivalent 
cytotoxic T lymphocyte (CTL) perforin-mediated killing, 
B6 mice are highly susceptible to PIFS, whereas 129S1 
mice remain resistant. This difference in susceptibility 
was illustrated by higher glial fibrillary acidic protein 
(GFAP) expression, increased CNS vascular permeability, 
and severe functional motor deficit in VP2i2i_i3o-treated 
B6 mice when compared to VP2i2i_i3o-treated 129S1 
mice. Furthermore, this susceptibility is transferable 
with the bone marrow compartment. 129S1 mice 
reconstituted with B6 bone marrow exhibit severe CNS 
vascular permeability, microhemorrhage formation, and 
functional deficit when compared to 129S1 mice receiv- 
ing autologous bone marrow transfer [14]. Through an 
initial microsatellite analysis, we have determined that 
the PIFS phenotype is not mediated through a single 
gene, but rather is likely to be a complex condition [15]. 

In order to provide a more comprehensive analysis 
and determine the location of genes that influence the 
PIFS phenotypes, we analyzed the F2 generation of 
brother-sister mated Fl animals from a B6 and 129S1 
cross using QTL analysis. The quantifiable traits of 
astrocyte activation via GFAP expression, CNS vascular 
permeability via fluorescent molecule diffusion across 
the BBB, and functional motor deficits measured with 
rotarod activity were chosen for study because they dif- 
fer by or exceed two standards of deviation when mea- 
sured in the two parental strains, which is optimal for 
QTL analysis [14,16]. We identified chromosomal re- 
gions that harbored QTLs responsible for functional def- 
icit and CNS vascular permeability. Additionally, we 
have identified candidate genes implicated in regulating 
these traits. These genes can be further investigated and 
may become potential targets for therapeutic approaches 
to ameliorate neurological disorders characterized by 
BBB disruption. 

Results 

As shown in Figure lA, B6 mice were crossed with 
129S1 mice to yield an Fl generation, and then Fl hy- 
brid mice were brother-sister mated to produce 303 F2 
progeny. These mice were intracranially infected with 
TMEV and then intravenously administered VP2i2i_i3o 
peptide 7 days later to induce PIFS. Mice were then 
analyzed for CNS vascular permeability, astrocyte 



expression of GFAP, and functional motor deficit (see 
Methods section). B6 mice were susceptible to each 
of these phenotypes, while 129S1 mice remained resist- 
ant [14]. F2 mice displayed variable susceptibility to 
each of these traits and were divided into susceptible 
and resistant groups based on whether they exhibited 
the trait or not. CNS vascular permeability was evaluated 
using a fluorescent plate reader to detect FITC-albumin 
leakage in brain homogenates using established methods 
[14, 17, IS]. Samples were read at A488, and a threshold 
value of 500 was used to determine susceptibility or re- 
sistance. Astrocyte expression of GFAP was detected 
through Western blot analysis as previously published 
[14, IS]. Expression values were based on GFAP protein 
levels normalized to GAPDH. We used a threshold value 
of 50 to categorize F2 mice into susceptible and resistant 
groups. Functional motor deficit was assessed using the 
rotarod behavioral assay. A threshold value of 50% initial 
motor ability was used to categorize mice as susceptible 
or resistant (Additional file 1). The phenotype data for 
each trait is illustrated by the mean ± SEM for the F2 
susceptible group and the F2 resistant group (Figure IB). 
No significant sex differences were observed in the F2 
mice (data not shown). 

Identification of QTL controlling quantitative traits 
associated with PIFS 

Of the 671 SNPs analyzed, the most significant linkage 
to susceptibility to PIFS was observed with markers on 
chromosomes 12 and 17, as illustrated by log-likelihood 
plots (Figure 2). Significant LOD scores were obtained 
for functional motor deficits, which mapped to chromo- 
some 12 (LOD = 3.3) (Figure 2A), and CNS vascular per- 
meability, which mapped to chromosome 17 (LOD = 3.7) 
(Figure 2B). No significant LOD score was obtained for 
astrocyte activation as measured by GFAP expression 
(Figure 2C). Parental allele distribution of the linked 
SNPs within the QTL intervals that mapped to the phe- 
notypes of functional motor deficits and CNS vascular 
permeability were detected using 2x3 contingency ta- 
bles comparing B6 homozygous (BB), heterozygous (BS), 
and 129S1 homozygous (SS) mice. The QTL location on 
chromosome 12 [PIFSl, LOD 3.3), estimated by a 95% 
confidence interval, was from 4.7 cM to 10.1 cM, with 
the peak of linkage at 6.0 cM. The QTL location on 
chromosome 17 {PIFS2, LOD 3.7), also estimated by a 
95% confidence interval, was from 5.7 cM to 12.5 cM, 
with the peak of linkage at lO.S cM. We found that the 
129S1 strain contributed to low rotarod performance, 
while the B6 strain contributed to high rotarod perform- 
ance on these markers (Table lA). In addition, the 
129S1 strain contributed to increased CNS vascular per- 
meability, whereas the B6 strain contributed to reduced 
CNS vascular permeability in the PIFS model (Table IB). 
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(A) 

C57BLy6J mice (PIFS -susceptible) 129S1/SvlmJ mice (PIFS -resistant) 



F1 (PIFS-resistant) F1 (PIFS-resistant) 



F2 generation 
(variable PIFS susceptibility) 

Analysis of: 1.) CNS vascular permeability 

2. ) Astrocyte expression of GFAP 

3. ) Functional motor deficit 



(B) 



Phenotype 


F2 Susceptible 


F2 Resistant 


CNS vascular permeability 
(detection of FITC -albumin 
leakage in brain homogenates) 


737.3 ± 23.2 

(^488) 


251 .7 ± 10.4 

lUppm (^488) 


GFAP expression 
(Western blot analysis) 


89.5 ± 3.2 
(GFAP:GAPDH ratio) 


22.3 ±1.4 
(GFAP:GAPDH ratio) 


Functional motor deficit 
(rotarod analysis) 


8.3% ± 1 .2% 
(% initial ability) 


157.5% ± 18.3% 
(% initial ability) 



Figure 1 Approach to map QTL for functional motor deficits, CNS vascular permeability, and astrocyte activation. (A) PIFS-susceptible B6 
mice were crossed with PIFS-resistant 129S1 mice to yield an Fl generation which was previously determined to be resistant to PIFS. Fl hybrid 
mice were then brother-sister mated to produce 303 F2 progeny with variable susceptibility to traits associated with PIFS, including functional 
motor deficits, astrocyte GFAP expression, and CNS vascular permeability measured by FITC-albumin leakage. DNA from 273 randomly selected F2 
progeny were analyzed using lllumina Chip Technologies. SNPs with a known cM value were analyzed for the relationship between quantifiable 
traits (functional motor deficits, CNS vascular permeability, and astrocyte expression of GFAP) and genotype. (B) F2 mice displayed variable 
susceptibility to each of the traits associated with PIFS and were divided into susceptible and resistant groups based on whether they exhibited 
the trait or not. CNS vascular permeability was evaluated using a fluorescent plate reader to detect FITC-albumin leakage in brain homogenates. 
Samples were read at A488, and a threshold value of 500 was used to determine susceptibility or resistance. Astrocyte expression of GFAP was 
detected through Western blot analysis. Expression values were based on GFAP protein levels normalized to GAPDH. We used a threshold value 
of 50 to categorize F2 mice into susceptible and resistant groups. Functional motor deficit was assessed using the rotarod behavioral assay. A 
threshold value of 50% initial motor ability was used to categorize mice as susceptible or resistant. All values are displayed as mean ± SEM. 
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Figure 2 QTL analysis of functional motor deficits, CNS vascular permeability, and astrocyte activation. Log-likelihood plots of (A) 
functional motor deficits as measured by rotarod performance, (B) CNS vascular permeability as measured by leakage of FITC-albumin into the 
CNS, and (C) astrocyte activation as measured by GFAP expression. The upper graph charts the data as a LOD score. The solid black line depicts 
the significance threshold. The lower graph depicts the estimated QTL effects. Significant LOD scores were obtained for (A) functional motor 
deficits, which is shown to map to chromosome 12 (LOD = 3.3), and (B) CNS vascular permeability, which is shown to map to chromosome 
17 (LOD = 3.7). 
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Table 1 Chi-square analysis of functional deficit and CNS vascular permeability 



A. SNP 




Low F2s (n = 


145) 






High F2s (n = 


128) 




Chi-square test 


Chi-square p-value 


LOD score 




BB 


BS 




SS 


BB 


BS 




SS 








rs6292954 


18 


80 




47 


50 


51 




27 


25.93 


0.000002 


3.65 


rsl 3481 303 


17 


83 




45 


46 


56 




26 


22.71 


0.000012 


3.66 


rs3655057 


19 


80 




46 


46 


54 




28 


19.66 


0.000054 


3.30 


rsl 3481 324 


21 


79 




45 


49 


52 




27 


20.28 


0.000039 


3.35 


B. SNP 




Low F2s (n = 


182) 






High F2s (n = 


:91) 




Chi-square test 


Chi-square p-value 


LOD Score 




BB 


BS 




SS 


BB 


BS 




SS 








rs6196216 


58 


88 




36 


20 


39 




32 


8.24 


0.016245 


3.73 


rs3672065 


60 


85 




37 


19 


39 




33 


9.27 


0.009706 


3.79 



2x3 contingency tables comparing B6 homozygous (BB), heterozygous (BS), and 129S1 homozygous (SS) mice with regards to (A) functional motor deficits as 
measured by rotarod performance and (B) CNS vascular permeability as measured by leakage of FITC-albumIn Into the CNS. (A) Mice displaying <50% Initial ability 
were placed In the low F2 group, while mice displaying >50% Initial ability were placed In the high F2 group. 129S1 contributes to low rotarod performance on 
chromosome 12 markers rs6292954, rsl 3481 303, rs3655057, and rsl 3481 324, while B6 contributes to high rotarod performance on these markers. (B) The low F2 
group consists of mice with FITC scores <500, while the high F2 group consists of mice with FITC scores >500. B6 contributes to less CNS vascular permeability on 
chromosome 17 markers rs6196216 and rs3672065, while 129S1 contributes to more CNS vascular permeability on these markers. 



Potential candidate genes implicated in regulating 
functional deficit and CNS vascular permeability in CDS T 
cell-initiated BBB disruption 

Using the Mouse Genome Informatics Database resources 
(www.informatics.jax.org), we performed a search between 
the base pairs of the significant loci found on chromosome 
12 (SNP markers rs6292954, rsl3481303, rs3655057, and 
rsl3481324) and chromosome 17 (SNP markers rs6196216 
and rs3672065) in order to discover genes located in prox- 
imity with these loci. Our search revealed 34 candidate 
genes on chromosome 12 and 58 candidate genes on 
chromosome 17 (Table 2). We then prioritized these genes 
based on known or presumed biological function and pre- 
vious studies, allowing us to narrow down this list to 6 
genes of interest on chromosome 12 and 8 genes of inter- 
est on chromosome 17. We focused on genes exhibiting 
potential immunological, hematopoietic, apoptotic, or 
CNS -related functions since factors related to these pro- 
cesses are likely to play a role in CD8 T cell- initiated BBB 
disruption. We then used the Mouse Phenome Database 
(www.phenome.jax.org) to determine the number of SNPs 
in these genes of interest that are different between the B6 
and 129S1 mouse strains as well as the location of these 
SNPs (Table 3). 

Discussion 

In this study, we investigated the role of genetic factors 
involved with CD8 T cell-initiated fatal BBB disruption. 
We previously demonstrated that despite having identi- 
cal MHC class I molecules, epitope-specific CNS- 
infiltrating CD8 T cells, and similar levels of perforin- 
mediated killing, the B6 strain is highly susceptible to 
PIFS, while the 129S1 strain remains resistant [14,15]. 
This susceptibility is transferable with the bone marrow 
compartment, as 129S1 mice reconstituted with B6 bone 
marrow exhibit characteristics of the PIFS phenotype. 



which include severe CNS vascular permeability, mic- 
rohemorrhage formation, and functional motor deficit, 
when compared to 129S1 mice receiving autologous 
bone marrow transfer [14]. Furthermore, by conducting 
a backcross of (129S1 x B6) Fl mice, we determined 
through microsatellite analysis that susceptibility to PIFS 
is dictated in a complex fashion [15]. Therefore, we in- 
vestigated the characteristics of PIFS in F2 generation 
mice through analysis of three quantifiable traits, which 
include functional motor deficit on rotarod, activation of 
astrocytes, and CNS vascular permeability. Since we 
have shown that the means of each of these traits mea- 
sured between the B6 and 129S1 mouse strains differ by 
greater than two standards of deviation, we expected to 
find a set of QTLs for each trait measured [14]. 

We identified one novel QTL (PIFSl) on chromosome 
12 linked to functional motor deficits and a second novel 
QTL {PIFS2) on chromosome 17 linked to CNS vascular 
permeability. While causality cannot be determined in this 
study, we have not observed motor dysfunction without 
CNS vascular permeability in the PIFS model. Further- 
more, we discovered that the SNP markers identified 
within the QTL intervals for the PIFS -resistant 129S1 
strain actually contributed to low rotarod performance 
and more CNS vascular permeability, while the B6 strain 
contributed to functional motor preservation and reduced 
CNS vascular permeability. While this was an unexpected 
result, it is possible that the expression of susceptibility 
loci may be suppressed by the presence of other loci in the 
parental strain. However, since new genetic combinations 
can be revealed through an F2 analysis, the identification 
of cryptic QTL becomes possible. For example, a genetic 
analysis of (B6 x FVB/NJ) F2 progeny to study atheroscler- 
osis revealed that the FVB/NJ strain, which is more re- 
sistant to atherosclerosis than the B6 strain, actually 
contributed to high levels of the disease [19]. Therefore, 
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Table 2 Candidate genes implicated in regulating Table 2 Candidate genes implicated in regulating 

functional deficit and CNS vascular permeability in fatal functional deficit and CNS vascular permeability in fatal 
BBB disruption BBB disruption (Continued) 



Gene 


Chr 


cM 


Genome coordinates 


Gene 


Chr 


cM 


Genome Coordinates 


NtSclb 


12 


5.27 


10376779-10396979 (+) 


Gml0513 


17 


7.8 


11918744-11925211 (-) 


Rdhl4 


12 


5.29 


10397586-10402368 (+) 


Gml6168 


17 


7.8 


10689312-10762262 (+) 


1700034J04Rik 


12 


5.58 


1 1228267-1 1229010 (-) 


Gml6169 


17 


7.8 


11000700-11019491 (+) 


493051 lA02Rik 


12 


5.58 


1 1445039-1 1451618 (+) 


Pacrg 


17 


7.8 


10595878-11033177 (-) 


Genl 


12 


5.58 


1 1247732-1 1272593 (-) 


Park2 


17 


7.8 


11033250-12256227 (+) 


Kcns3 


12 


5.58 


1 1097008-1 1 157648 (-) 


Agpat4 


17 


8.33 


12311570-12412511 (+) 


Msgnl 


12 


5.58 


1 1215188-1 1215754 (-) 


Map3k4 


17 


8.42 


12420487-12511526 (-) 


Rad51ap2 


12 


5.58 


1 1462885-1 1469734 (+) 


4732491 K20Rik 


17 


8.5 


12511751-12520116 (+) 


Smc6 


12 


5.58 


11272692-11326591 (+) 


Pig 


17 


8.5 


12571475-12612250 (+) 


Vsnll 


12 


5.58 


11332053-11443455 (-) 


Slc22a3 


17 


8.52 


12612838-12700570 (-) 


Fam49a 


12 


5.63 


12268945-12383165 (+) 


C030013G03Rik 


17 


8.56 


12659710-12672444 (+) 


Gml7337 


12 


6 


12794791-12801091 (+) 


Slc22a2 


17 


8.61 


12777055-12821331 (+) 


493051 9A1 IRik 


12 


6.1 1 


12909526-12914057 (+) 


Slc22al 


17 


8.63 


12841740-12868704 (-) 


Mycn 


12 


6.14 


12942902-12948720 (-) 


Igf2r 


17 


8.64 


12875272-12962530 (-) 


Ddxl 


12 


6.4 


13226113-13255980 (-) 


Aim 


17 


8.66 


12934177-13052988 (+) 


Nbas 


12 


6.45 


13275933-13590616 (+) 


Gml9270 


17 


8.66 


12926941-12946579 (+) 


Fam84a 


12 


6.93 


14154405-14158844 (-) 


Gml7592 


17 


8.69 


13030801-13031000 (-) 


Gm 16497 


12 


7.03 


14268340-14292490 (+) 


Masi 


17 


8.69 


13033945-13061009 (-) 


Trib2 


12 


7.3 


15798533-15823683 (-) 


Mrgprh 


17 


8.7 


13068900-13070708 (+) 


Gm5432 


12 


7.28 


15823809-15894961 (+) 


Pnldcl 


17 


8.71 


13081768-13102866 (-) 


Lpin1 


12 


7.9 


16542475-16596576 (-) 


Acat3 


17 


8.72 


13116699-13133268 (-) 


Ntsr2 


12 


7.98 


16660276-16667042 (+) 


Mrpll8 


17 


8.72 


13104215-13109211 (-) 


Grebl 


12 


8 


16677421-16807692 (-) 


Snora20 


17 


8.72 


13115723-13115783 (+) 


2410004P03Rik 


12 


8.04 


1701 1764-17018360 (-) 


Tcpl 


17 


8.72 


13108567-13117933 (+) 


E2f6 


12 


8.04 


16817771-16833549 (+) 


Acat2 


17 


8.73 


13135756-13153613 (-) 


Kcnfl 


12 


8.04 


17178906-17183694 (-) 


Wtap 


17 


8.73 


13159662-13185412 (-) 


Pqlc3 


12 


8.04 


16995454-17007214 (-) 


Sod2 


17 


8.75 


13200705-13210985 (+) 


Rock2 


12 


8.04 


16901784-16995080 (+) 


Gpr31b 


17 


8.76 


13244187-13245146 (-) 


Atp6vlc2 


12 


8.06 


17291527-17336165 (-) 


Tcpl Ob 


17 


8.76 


13253970-13275347 (+) 


Pdia6 


12 


8.06 


17273351-17291576 (+) 


Gmlll66 


17 


8.78 


13289526-13291338 (-) 


NollO 


12 


8.07 


17355299-17436901 (+) 


Unc93a 


17 


8.78 


13301483-13324657 (-) 


Mir3066 


12 


8.08 


17362198-17362280 (+) 


Gml0512 


17 


8.81 


13397908-13399077 (+) 


Gml9657 


12 


8.09 


17436584-17441442 (-) 


Smok2a 


17 


8.82 


13414054-13420524 (+) 


Odd 


12 


8.11 


17551679-17558308 (+) 


Smok2b 


17 


8.82 


13421718-13430055 (+) 


PdelOa 


17 


4.96 


8718237-9179513 (+) 


TcplOc 


17 


8.86 


13547438-13570089 (+) 


170001 01 14Rik 


17 


5.53 


9181198-9201184 (+) 


Rnu6 


17 


8.88 


13624899-13625005 (-) 


653041 IMOlRik 


17 


5.97 


9340584-9360701 (-) 


2700054A10Rik 


17 


8.9 


13675419-13898832 (-) 


Gml7728 


17 


7.33 


9614925-9615323 (+) 


Tcte2 


17 


8.9 


13675419-13898832 (-) 


Gm 17748 


17 


7.33 


9614805-9622259 (-) 


Gml6050 


17 


8.91 


13695065-13704059 (+) 


Pabpc6 


17 


7.47 


9859362-9862569 (-) 


Gm8603 


17 


8.91 


13709618-13711010 (+) 


Qk 


17 


7.75 


10399336-10512226 (-) 


Smok4a 


17 


8.91 


13714322-13721299 (+) 


A230009B12Rik 


17 


7.8 


10649082-10685921 (+) 


4930488N24Rik 


17 


8.95 


14238826-14243457 (-) 



Johnson et at. BMC Genomics 2013, 14:678 
http://www.bionnedcentral.conn/1471 -21 64/1 4/678 



Page 7 of 10 



Table 2 Candidate genes implicated in regulating 
functional deficit and CNS vascular permeability in fatal 
BBB disruption (Continued) 



Dact2 


17 


8.95 


14332238-14340838 {-] 


Gm 16046 


17 


8.95 


13812667-13820523 {-] 


Gm 16049 


17 


8.95 


13896413-13897375 {-] 


Gml6052 


17 


8.95 


13896377-13896992 {-] 


Gm7168 


17 


8.95 


14085380-14087685 (+ 


Gm7356 


17 


8.95 


14137246-14138709 {-] 


Gm7358 


17 


8.95 


14195012-14196538 {-] 


Mllt4 


17 


8.95 


13897546-14043157 



+) 



Shown are a list of genes between the significant loci found on chromosome 
12 (SNP markers rs6292954, rsl 3481 303, rs3655057, and rsl 3481 324) and 
chromosome 17 (SNP markers rs6196216 and rs3672065). Genes that 
exhibited potential immunological, hematopoietic, apoptotic, or 
CNS-related functions were categorized as genes of interest. 

while our finding may seem surprising, it is important to 
note that disease-resistant strains also harbor susceptibility 
alleles, consistent with other genetic studies [19-21]. 

Identification of QTL regions that influence resistance 
and susceptibility to PIFS has enabled us to identify candi- 
date genes that mediate protection from PIFS and thus are 
potential targets for therapy in diseases characterized by 
BBB disruption. Genes whose biological function had rele- 
vance to playing a role in BBB disruption were categorized 



as genes of interest and will be targeted first in future stud- 
ies. Genes of interest pertaining to mouse functional ability 
on chromosome 12 include tribbles homolog 2 {Trib2), 
which has been shown to induce apoptosis of cells in 
hematopoietic origin. Additionally, mice homozygous for a 
null Trib2 allele display increased mean percentage of 
CDS T cells. Furthermore, Trib2 has also been character- 
ized as one of the top non-HLA markers associated with 
susceptibility to multiple sclerosis [22]. Considering our 
previous research on the critical roles of CDS T cells and 
hematopoietic factors in mediating BBB disruption, Trib2 
will be a candidate for future study [9,10,14,15,17,1S]. 
Other genes of interest include Lpinl and Rock2, which 
both play a role in organization of the actin cytoskeleton, 
Vsnll, which plays a role in neuronal Ca^^ sensor proteins, 
Kcnfly which plays a role in neurotransmitter release and 
neuronal excitability, and Ddxly which plays a role in viral 
responses. Genes of interest on chromosome 17 that may 
affect CNS vascular permeability include PdelOa, Q/c, 
Pacrg, Park2, Pig, Igf2n Masl, and Mllt4, These genes ex- 
hibit a variety of functions that may be involved in causing 
disruption of the BBB and ensuing CNS vascular perme- 
ability. For example, Mllt4 has been shown to play a role 
in adherens junctions and cell junctions. Pig has been 
demonstrated to be involved in apoptotic processes and 
vessel development. Additionally, Masl is involved in 



Table 3 Genes of interest potentially related to B6 and 129S1 phenotypic differences in CDS T cell-initiated BBB 
disruption 



Chromosome 12 (Rotarod) 


Proposed function 


# SNPs different between 
strains 


SNP locations 


Vsnll 


Neuronal Ca2+ sensor proteins 


24 


All 1 


Ddxl 


Viral responses 


2 


All 1 


Trib2 


Induces apoptosis of cells in hematopoietic origin; 
CDS T cells 


40 


1 Cs, 1 U5, 38 1 


Lpinl 


Demyelination; actin cytoskeleton reorganization 


531 


3 Cn, 8 Cs, 27 U3, 
493 1 


Kcnfl 


Neurotransmitter release; neuronal excitability 


35 


4 Cs, 7 U5, 24 U3 


Rock2 


Actin cytoskeleton organization; regulation of 
angiogenesis 


330 


1 Cn, 3 Cs, 5 U3, 
321 1 


Chromosome 17 (CNS vascular 
permeability) 


Proposed function 


# SNPs different between 
strains 


SNP locations 


PdelOa 


Signal transduction; neuronal cell bodies 


1268 


Unspecified 


Qk 


Myelinization; axon ensheathment; vasculogenesis 


214 


12 U3, 202 1 


Pacrg 


Myelinization 


265 


2 U5, 263 1 


Park2 


Neuron projection; neuron death 


3198 


2 Cs, 1 U3, 3195 1 


Pig 


Apoptotic processes; vessel development 


38 


1 Cn, 37 1 


Igf2r 


Edema 


7 


All 1 


Masl 


Inflammatory responses 


2 


All 1 


Mllt4 


Adherens junctions; cell junctions 


7 


All 1 



Shown are the genes of interest found on chromosome 12 and chromosome 17 along with their proposed functions. The Mouse Phenome Database 
(www.phenome.jax.org) was used to determine the number of SNPs that are different between the B6 and 129S1 mouse strains as well as their location. 
{Abbreviations: Cn nonsynonomous codon, Cs synonomous codon, U5 5^ UTR variant, U3 3^ UTR variant, / Intron-varlant). 
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inflammatory responses, IgfZr is involved in edema, and 
Qk and Pacrg both play a role in myelinization. Further- 
more, Park2 has been shown to contribute to neuron pro- 
jection as well as neuron death, and PdelOa has been 
shown to be involved in signal transduction and neuronal 
cell bodies (www.informatics.jax.org/). Furthermore, we 
used the Online Mendelian Inheritance in Man (OMIM) 
database (http://omim.org/) to investigate whether any of 
these genes of interest have been associated with any 
human neurological diseases. We found that Vsnlly 
which is located on mouse chromosome 12, has been 
implicated in Alzheimer's disease. Q/c, located on mouse 
chromosome 17, has been proposed to play a role in 
myelin and oligodendrocyte dysfunction in schizophre- 
nia. Additionally, Pacrg, also located on mouse chromo- 
some 17, may be involved in acquired, communicating 
hydrocephalus (http://omim.org/). 

Our comprehensive genetic analysis has allowed us to 
create this list of candidate genes, which will serve as a 
critical starting point in elucidating the underlying mo- 
lecular mechanisms of BBB disruption in neurological 
diseases. This will be an essential first step in developing 
potential therapeutic approaches to treat these diseases. 
Validating strain differences in primary DNA sequences 
in candidate genes and performing bioassays to test for 
functional consequences of primary DNA differences 
will be essential in this process. It will also be essential 
to test candidate genes for their ability to contribute to 
susceptibility or resistance to functional deficit and BBB 
disruption in the PIFS model. Testing specific congenic 
strains in which potential genes of interest will be intro- 
duced will help define the importance of the QTLs 
detected in this study. This approach has been successful 
in other genetic studies [23,24]. Generation of bacterial 
artificial chromosome (BAG) transgenic mice using the 
candidate genes we have identified will also be beneficial 
in evaluating the importance of our novel QTLs. This 
will enable us to test whether strain-specific variation in 
one or more of these genes is responsible for causing the 
QTL effects we observed [25]. This method has been 
carried out successfully in other genetic studies of strain 
comparison, such as evaluating the difference in seizure 
susceptibility between B6 and DBA/2 mouse strains [25]. 
Ultimately, results obtained from these experiments may 
hold promise for therapeutic intervention of immune- 
mediated neurological diseases characterized by exten- 
sive BBB disruption. 

Conclusions 

The comprehensive genetic analysis performed in this 
study has enabled us to identify QTL linked to func- 
tional motor deficits and GNS vascular permeability, 
which are both devastating characteristics associated 
with GD8 T cell-initiated BBB disruption. Furthermore, 



identification of these critical QTL intervals has revealed 
a list of candidate genes implicated in regulating these 
traits. These findings are important because they will set 
the stage for future research to identify novel factors 
encoded by these genes that affect immune-mediated 
BBB disruption. This may ultimately hold promise for 
devising potential therapeutic approaches to treat the 
numerous immune-mediated neurological disorders 
characterized by BBB disruption. 

Methods 

Animals 

Male and female G57BL/6J (B6) (stock# 000664) and 
129Sl/SvImJ (129S1) (stock# 002448) mice were 
obtained from Jackson laboratories (Bar Harbor, ME). 
(B6 X 129S1) Fl hybrids were brother-sister mated to 
produce (B6 x 129S1) x (B6 x 129S1) F2 progeny. All Fl 
and F2 mice were produced in the animal colony at the 
University of Gincinnati (Gincinnati, OH). In order to 
obtain large cohorts, mice were tested at 4-7 months of 
age. All experiments were approved by the Institutional 
Animal Gare and Use Gommittee of the University of 
Gincinnati and Mayo Glinic. 

Induction of CNS vascular permeability using the PIFS 
model 

GNS vascular permeability was induced as previously de- 
scribed [15,17,18]. Briefly, all mice were intracranially 
infected with 2 x 10^ PFU Daniels strain of TMEV. 7 days 
post-TMEV infection, during the peak of GD8 T cell ex- 
pansion, mice were intravenously administered 0.1 mg 
VP2i2i.i3o (FHAGSLLVFM) peptide [15]. 

Rota rod 

In order to evaluate functional motor deficits, mice were 
placed on the Rotamex-5 rotarod apparatus (Golumbus 
Instruments) increasing from 5-40 RPMs over 7 minutes. 
Mice were trained twice daily for 3 consecutive days. The 
motor performance of mice was then assessed the day fol- 
lowing this training period in order to establish baseline 
ability, which was considered the average of two trials. 
Mice were then injected with VP2i2i-i3o peptide to induce 
GNS vascular permeability and their motor performance 
was assessed again 24 hours later. Final scores were 
depicted as percent initial motor ability (day 2 average/day 
1 average x 100). 

Fluorescein isothiocyanate (FITC)-albumin permeability assay 

Mice were given an intravenous injection of 10 mg 
FITG-albumin (Sigma #A9771) 23 hours post-VP2i2i_i3o 
peptide administration to induce PIFS. Whole brains 
were harvested one hour after administration of FITG- 
albumin and frozen on dry ice. Brain tissue samples were 
homogenized in radioimmunoprecipitation assay (RIPA) 
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buffer [10 mmol/L Tris, 140 mmol/L NaCl, 1% Triton 
X-100, 1% Na deoxycholate, 0.1% SDS and protease in- 
hibitor cocktail (Pierce #78410) pH 7.5] and spun for 
15 min at 10,000 RPM. Protein concentration was 
assessed using the BCA protein assay (Pierce #23223). 
Samples were then normalized for protein content and 
read on a fluorescent plate reader at 488 nm excitation 
and 525 nm emission to detect FITC-albumin leakage 
into the brain. Data was collected using SpectraMax 
software (Molecular Devices). 

Western blot for glial fibrillary acidic protein (GFAP) 
detection 

Normalized protein samples from the FITC-albumin per- 
meability assay were also used for western blot analysis. In 
order to detect GFAP, 10 ug of protein was loaded per well 
on 4-20% precise protein gels (Pierce #25244) with BupH 
Tris-HEPES-SDS running buffer (Pierce #28398). Tris 
transfer buffer [400 mmol/mL Tris base, 70 mmol/mL gly- 
cine, 10% methanol] was used to transfer gels onto 
Immun-Blot PVDF membranes (Biorad #162-0177). GFAP 
was detected by staining with mouse anti-GFAP (1:1000 
BD Pharmingen #556329) and goat anti-mouse IgG conju- 
gated to horseradish peroxidase (Sigma #A3682). Scion 
Image Software (Scion Corp) was used to analyze Western 
blot films. Values were based on GFAP protein levels nor- 
malized to GAPDH [14,18]. 

Quantitative trait linkage analysis 

DNA from 273 F2 progeny (124 males; 149 females) 
were analyzed on 2 SNP chips using Illumina Chip 
Technologies. Two B6 and 2 129S1 mice were included 
as controls. In order to map genes controlling suscepti- 
bility to PIFS in the B6 x 129S1 F2 model, we 
constructed a genetic linkage map of the mouse genome 
by analyzing 273 randomly selected mice using Illumina 
Chip Technologies. Of the 1449 analyzable SNPs, 874 
were found to be unique between the two strains. 671 of 
these SNPs had a known cM value and thus were ana- 
lyzed for the relationship between quantifiable traits and 
genotype [26]. Based on previous studies, we used func- 
tional motor deficits measured by rotarod, CNS vascular 
permeability measured by FITC-albumin leakage, and 
astrocyte activation measured by GFAP expression as 
quantifiable traits predictive of susceptibility to PIFS 
[14,15]. Threshold values were generated using the 
WinQTL Cartographer with 1,000 permutations of our 
data set to correct for multiple testing [27]. Interval 
mapping was used to identify significant linkage using a 
critical value of p = 0.05 and a 95% confidence interval 
(Additional file 2). A chi-square test statistic for each 
SNP marker showing linkage within a QTL interval was 
obtained using 2x3 contingency tables to demonstrate 
the distribution of parental strain alleles within the F2 



population using the quantitative traits of functional 
motor deficits, CNS vascular permeability, and astrocyte 
activation. Since F2 mice displayed variable susceptibil- 
ity, they were divided into two groups, low F2s and high 
F2s, based on whether or not they exhibited functional 
motor deficits or CNS vascular permeability. For func- 
tional motor deficits, low F2s were defined as those mice 
displaying <50% initial ability, while high F2s were de- 
fined as mice displaying >50% initial ability. For CNS 
vascular permeability, mice with FITC scores <500 were 
placed in the low F2 group while mice with FITC scores 
>500 were placed in the high F2 group for analysis. Since 
data were scored in a binomial fashion, we used a chi- 
square test, which is a nonparametric approach, to de- 
termine significance. We identified the critical interval 
for significant QTL and performed a search using www. 
informatics.jax.org to compile a list of candidate genes. 
These genes were then prioritized based on their known 
or presumed biological function. We then used the 
Mouse Phenome Database (www.phenome.jax.org) to 
determine the number of SNPs in these genes of interest 
that are different between the B6 and 129S1 mouse 
strains as well as the location of these SNPs. These genes 
will set the stage for future experiments on gene expres- 
sion and strain differences that may be related to the 
mechanism of CD8 T cell-initiated BBB disruption. 

Availability of supporting data 

The data sets supporting the results of this article are in- 
cluded as additional files. 

Additional files 



Additional file 1: Shown are the raw data obtained for CNS 
vascular permeability, as measured by assaying FITC-albumin 
leakage in brain homogenates, astrocyte expression of GFAP, as 
detected by Western blot analysis, and functional motor deficit, as 
assessed using the rotarod behavioral assay. Mice were divided into 
susceptible and resistant groups using the criteria stated in the text. Raw 
data values were used to determine mean ± SEM for each trait. 

Additional file 2: Shown are the raw data used to analyze the SNPs 
with known cM values for the relationship between quantifiable 
traits and genotype using WinQTL Cartographer. Mouse identification 
numbers are shown across the top. SNP markers, chromosome numbers, 
cM values, and genotype data are included. For genotype, 0 = 
homozygous for B6, 1 = heterozygous, and 2 = homozygous for 129S1. 
Raw data for each quantifiable trait are included in a separate tab. 
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